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Abstract 

A general input-output modelling technique for aperiodic-sampling lin- 
ear systems has been developed. The procedure describes the dynamics 
of the system and includes the sequence of sampling periods among the 
variables to be handled. Some restrictive conditions on the sampling 
sequence are imposed in order to guarantee the vahdity of the model. 
The particularization to the periodic case represents an alternative to the 
classic methods of discretization of continuous systems without using the 
Z-transform. This kind of representation can be used largely for identifi- 
cation and control purposes. 



1 Introduction 

Aperiodic sampling is a very interesting technique for improving the sohition of 
several problems in control and identification. In particular, aperiodic sampling 
systems have been found useful in signal adaptation and compensation [4l [2] and 
in optimal transmission of measuring errors in problems involving the solutions 
of systems of linear equations, such as observability [HIS], controllability [8] and 
identifiability [7l|8]. 

In all these cases it is essential to find a general model that 

(a) describes the dynamics of the system; 

(b) is adapted to real experimentation conditions (external representation against 

internal-state representation); 

(c) includes the sequence of sampling periods among the variables to be han- 

dled. 
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A well-known input-output formulation is used successfully for linear systems 
sampled periodically. The scalar output of the plant at an arbitrary instant is 
described by means of input and output samples taken at previous instants. It 
seems natural to consider the aperiodic case in the same way, except that the 
finite difference equation coefficients, depending on the fixed sampling period, 
would have to be replaced by multivariable functions, depending on the aperi- 
odic sampling sequence. Such a formulation would satisfy the above conditions. 

There must be restrictions on the aperiodic sampling sequence in order that 
the model be valid, whereas in the periodic case any sampling interval is valid. 

The systems considered are characterized by their impulse response (weight- 
ing function), which describes the output of the plant through a convolution 
expression. Consequently, the properties and results obtained are always given 
in terms of that function. 

The paper is organized as follows. In §1 several basic assumptions and gen- 
eral considerations, which will be assumed through the work, are given. Section 
2 is devoted to a modelling technique for aperiodic-sampling linear systems. 
This also includes the action of a sampler and zeroth-order hold preceding the 
plant. Section 3 contains some restrictive conditions on the sampling sequence. 
The problem of the selection of the aperiodic sequence is strictly considered in 
geometric terms and the results found in the literature are improved. In §4 a 
particularization to the periodic case has been made. The proposed technique 
represents an alternative to the classical methods of discretization of continuous 
systems, in order to obtain the discrete-model coefficients, without using the 
Z-transform. Finally conclusions in §5 end the paper. 

2 Basic assumptions 

This discussion is restricted to: 

(i) linear time- invariant single- input / single-output differential systems of finite 

order n; 

(ii) systems whose transfer fimction G{S) is a strictly proper rational function, 

and whose impulse response h{t) is therefore a particular solution of 

/i(")(t) + ai/i("-i)(i) + ... + a„/i(t) = 0, <>0 (1) 

an nth-order homogeneous linear differential equation with constant coef- 
ficients {tti e M). 

h{t) can be then be written as 

n 

h{t) = j2a(p,{t), t>o (2) 

where Ci € C are constant coefficients and (^i; M — > C (z = 1, ri) is the 
fundamental system of solutions of eqn. ([1]). 
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We conclude this preliminary section with the following statement. Let (Gi) 
be a family of vector functions 

Gi:M"^M" (i = 0, !,...,«) 

Gi e C°°(M",M") 

C°°(R", R") being the set of infinitely differentiable functions on R". It then 
follows that if there exist 

(a) an integer r < n such that the elements {Go{z), . . . ,Gr{z)) are linearly 
independent for all z e R", and 

(6) an integer k > r such that Gk{z) depends linearly on {Go{z), . . . ,Gr{z)), 
then there are functions 

/o,/i,...,/neC°°(M",M) 

such that 

n 

^/,(z)G,(z) = VzeR" (3) 



3 Modelling technique 

3.1 Generalities and key ideas of the methodology 

The class of linear time- invariant SISO systems is characterized by the impulse 
response, which describes for zero initial conditions the output of the plant 
through a convolution expression 



Vk 



(4) 



i=0 



where yk is the output of the plant at time tk, Ui is the impulse input at time 
U, and h{t) is the impulse response of the plant. According to ([1]), the impulse 
response can also be written in matrix form by means of the equivalent linear 
system. 



X{t) = AX{t) 



(5) 



where 



xit)^[h{t),hit),...,h(--'Ht)r (6) 

the symbol ' denotes the transpose and A is an n x n bottom-companion matrix. 



A 



(7) 



-a„ — a„_i —0,1-2 



-ai 
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The solution of the linear system 

X{t) = exp{At)Xo {Xo = X{0)) (8) 
is related to the impulse response through the expression 

h{t) = c exp{At)Xo (9) 

with 

c=[l,0,...,0] (10) 
Note that the components of the vector Xq 

Xo = [huh2,...,hn]' (11) 

correspond to the first n Markov parameters 

d'h{t) 



h, 



+1 



(i = 0,l,...,n-l) (12) 

t=o 



It should also be noted that the triad {A,Xo,c) leads us naturally to the 
observability canonical realization from the scalar impulse response. 

From these matrices A,Xq,c we are going to define a family of vectors func- 
tions 

(Go, Gi, . . . , Gn) 
Gi:M"^-M" (i = 0,l,...,n) 

given by 

Gi{z-i,..., Zn) = exp{A{zi + . . . + Zi))Xo {i = l,...,n) (13) 



Go{zi,...,Zn)=Xo (14) 

with 

The functions Gi and the impulse response are related by means of the 
expression 

h{zi + ... + Zi) = cGi{z) zi + ... + Zi>0 (15) 

Let us now consider the statement given in §1 for the kind of function defined 
above. 

Prom an analytical viewpoint, the functions Gi belong to C'°°(M",M"), as 
compositions of C°° functions. 

Let /" C K" be an open subset of M" such that the vectors {Go{z),... 
, G„_i(z)) are linearly independent for all z G In this case, it is easy to see 
that for the new domain 7" the conditions (a) and (b) in the previous statement 
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hold. In fact, condition (a) holds by definition of the subset and condition 
(h) holds by the dimensionality of the vector Gi. Hence there will be functions 



e C°°(/",M) (z = 0,...,n) 



such that 



Y,fn-^{z)G,{z)^Q Vze/" 

Rewriting (jl6p in matrix form, with 

G^{z) = [Gu{z) . . . Gm{z) ]' = 0, . . . , n - 1) 

we get 



(16) 



Gio{z) . 


.. Gi,„_i(z) ■ 




' Uz) ' 




" Gi„(z) " 


Gno{z) ■ 


Gn,n-l{z) 








Gnn{z) 



(-/oW) (17) 



{fo{z) = —1 for simplicity) and the functions (/i(z), . . . , fn{z)) can be ob- 
tained by solving a compatible system of linear equations. 
The general form of the functions fi{z) is 



Mz) 



Det[Go{z)...Gn{z)...Gn-i{z)\ 



(18) 



Det[Go{z)...Gn-i{z)] 
where the numerator is the determinant obtained from the matrix 

[Go(z),...,G„_i(z)] 

by replacing the zth column by the column vector G„(z). 

At this point, we identify the components of z with the elements of the 
sampling period sequence 



Zn =tk — tk-1 — Tk 
Zn-1 — tk-1 — tk-2 — Tk-1 

Zl — tk-n+1 — tk-n = Tk-n+1 



(19) 



We multiply both sides of ([TE)) by cexp{Az*), with z* successively taking 
the values 



■^zi (i = l,...,n) 



(20) 



1=1 
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and we get in each case 

n — i 

cexp{A{zi+i + . . . + Zn))Xo = c J2 fi{z) exp{A{zi+i + . . . + z„_;))Xo+ 

1=1 

i-l 

c J2 fn-i{z) exp{~A{zi+i + ... + Zi))XQ 

{i = l,...,n) 

(21) 

We define 

gn^i{z) = c^fn-i{z)exp{A{zi+i + ... + z^))Xo {i = l,...,n) (22) 

1=0 

From ((22|) the functions gn-i{z) are of class C°° as compositions of C°° 
fmictions. 

It can be checked that the functions gn~i{z) correspond to the first compo- 
nent of hnear combinations of vectors Gi with negative argument. 
If we write 

= m 

9j ^djiz) 

at time tk , we can condense the preceding expressions into two sets of equa- 
tions involving the functions fi,gj and h{t): 

n 

Y^J^h{tk-^-tj)+gl_^={) (j-fc,fc-l,...,fc-n + l) (23) 

n 

^ /f /i(ife_, - tj) =0 (fc > n; J = fc - n, . . . ,0) (24) 

i=0 

Multiplying every equation by Uj (j = fc, fc — 1, . . . , 1, 0) (impulse inputs at 
the sampling instants) respectively and summing, we get 

k—l k—n 
f'l ^ Ui h{tk^i ~ti) + ...+f^^Ul h{tk-n - tl) 
1=0 1=0 

n — 1 k 

+ ^g''jUk-j ^^uih{tk-ti) (25) 

3=0 1=0 

and, according to Q, the preceding expression becomes 

n n—1 

Vk^Y. ^^y^-' + H 9jUk-j (26) 
i=l j=o 

which is called the input-output model for linear time-invariant aperiodic- 
sampling systems. 
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3.2 Simplified form of the functions fi,gj 

Companion matrices are an important example of what are called cyclic (or non- 
derogatory) matrices, which have only one (normalized) eigenvector associated 
with each distinct eigenvalue. This means that 

(i) the Jordan canonical form is clearly simplified (there is only one Jordan 

block for each distinct eigenvalue); 

(ii) the similarity transformation of the given matrix A to the Jordan canonical 

form can be obtained in a standard way. 

Indeed, 

A = BJB^^ (27) 

where J is the Jordan canonical form of the matrix A, and B is an invertible 
matrix of a well-known general form. 
In this way, (|16p becomes 

n 

J2fn-r{z)exp{Ja,)Yo = (28) 



i=0 



with 



a, = zi + . . . + Zi (i = 1, . . . , n) (29) 
ao = (30) 

Yo = B-'Xo (31) 

Factorizing det[exp{JaQ)Yo, . . . , ( Ja„_i)yo] by means of Laplace's expansion 
by minors and cancelling common factors in the numerator and denominator of 
(fTB]), the functions can be written as 

ft = ^ (z = l,...,n) (32) 

where 

A^\ipi{a^)\ (i = 0,l,...,n-l) ] 
A,^\^i{a,)\ (j = 0,...,n,...,n-l) \ (33) 
(1^1,. ..,n) J 

The linear independence of the vectors Gj implies the non-nullity of the 
determinant A. We can check that the functions /j depend exclusively on the 
poles of the transfer function and on the sampling sequence. The result could 
be expected and agrees for the periodic case with the direct correspondence 

Xi -J> exp{XiTo) 

between the poles of the continuous and pulse transfer function. 
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The functions g'^ are given by 



9', 



-j=-Y.^^^^^k-^-^i) ij =k,k~l,...,k-n + l) (34) 



and depend on the poles and zeros of the transfer function, since they are 
obtained as a hnear combination of the impulse response for specific arguments. 

3.3 General formulation for systems with zeroth-order hold 

If the sampler at the input is followed by a zeroth-order hold the input-output 
model is modified in the following way: 

k 

yk = ^h{tk-ti)xi (35) 

1=0 

with 

Xi =Ui - iti_i (36) 

h{t) = £-1 [ 1 (37) 

As Xi represents the difference between two consecutive values of the input 
signal, the expressions (|23l) and (fM)) become 

n ^ 

E f^[hitk-^ - t,) - h{tk-, - i,+i) ] + gt, =0 {j ^k,k~l,...,k-n) 
£ ff[h{tk^, - t,) - h{tk-^ - tj+i) ] = (j = fc - n - 1, . . . , 0) 

i=0 

(fc > n) 

(38) 

Briefly, the presence of a zeroth-order hold implies that 
(a) the functions /,; are the same as before; 

(&) the functions gj are modified: the number functions gj is increased by one 
(so that now < j < n), and their general expression is similar to (|34p 
but with 

h{tk-i — tj) h{tk_i — tj) — h{tk-i — tj+i) 
for h{.) defined in ([57|). 

The results obtained were to be expected, since the functions fi depend 
only on the poles of the transfer function, which reflect the internal coupling 
in the system and its autonomous behaviour. However, the functions gj reflect 
the internal plant coupling to the input signal, which has been affected by the 
presence of the zeroth-order hold. 
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3.4 Main results 

We recall briefly the main results we have obtained. 

(a) We have shown that the functions gj are infinitely differentiable. 

(6) We have obtained a general and systematic formulation for every function 
fi, gj in contrast with the procedure found in [6 , where all these functions 
are computed globally. 

(c) We have developed a procedure to impose restrictive conditions on the 
sampling sequence in order to guarantee the linear independence of the 
vectors (Go(z), Gi{z), Gn-i{z)). 

In the next section we shall determine the set of vectors /" C M" whose elements 
z satisfy the above condition. 

4 Choice of the samphng-period sequence. A 
geometric interpretation 

The problem of the choice of the sampling sequence has been treated in [9] 
analytically. Making use of the concept of a Chebyshev system, some intervals of 
the real line are selected in which the sampling instants can be chosen freely. In 
this way, the non-nullity of the determinant A is guaranteed, and consequently 
the difference equation for the aperiodic case can be obtained (via (1^ . ((5^ 



In the present paper the same problem is considered geometrically. The 
choice of sampling-period sequence is directly related to the properties of certain 
vectors in the space M". In this way, the intervals found by Troch [5] are largely 
increased for low-order models and some general considerations for higher-order 
models are made. 

4.1 Second-order model (n = 2) 

We are going to consider a second-order model with a pair of complex eigenvalues 



The problem depends on an adequate choice of the sampling periods in such 
a way that the linear independence of the vectors 



and (p4|). 



a + jbeC 



(6>0) 



(39) 



[exp{JaQ) zq, exp{Jai) zq] 



(40) 



will be preserved, where 



ai = tk-i — tk-2 — Tk-i 



(41) 
(42) 



Jr is the real canonical form of the matrix A: 
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(b > 0) (43) 
and T is the (invertible) matrix governing the change of basis: 

zo = T-^Xo (44) 

As we are in M^, the geometric interpretation is very simple. The generic 
operator exp{Jra) apphed to the vector zq can be viewed as foUows. It is 
a counterclockwise rotation through ba radians, followed by a stretching (or 
shrinking) of the length of zq by a factor exp{aa) [3]. 

From this interpretation Yq and Yi, will be linearly independent if and only 

if 

bai^bTk-i^n (45) 

where n denotes an integral multiple of vr. Otherwise the vectors will be colinear. 

Comparing these results to those of Troch j9] we have the following. Ac- 
cording to [9], given ifc-2, ifc-i can be fixed so that 

tfc-i e (ife-2,ifc-2 + ^) (46) 
According to (|45p. given tk-2, ife-i can be fixed so that 

tk-i e / = (ifc-2, oo) - (tk-2 + ^) (47) 
Only point values will be rejected. 

It can also be checked that such a formulation is concerned with differences 
between sampling instants but not with the actual position of such instants on 
the real axis. This agrees perfectly with the kind of time-invariant systems we 
are dealing with. 

If the condition pS)) is violated, the sampling process resonates with the 
system dynamics and the formulation then obtained no longer affords a faithful 
representation of the system. 

We can also see that, for the second-order model, (j45|) is a necessary and suf- 
ficient condition to guarantee the linear independence of the vectors Yi whereas 
the condition given by Troch is only sufficient. 

4.2 Third-order model (n = 3) 

We are going to consider a third-order model with a real pole and a complex 
pair 

A G M, a + jbeC {b>0) (48) 

The point is to choose the sampling instants in the right way; that is, such 
that the vectors 



Jr = 
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[^0,^1,^2] = [exp{Jrao) zo,exp{Jrai) zo,exp{Jra2) zq] (49) 
are linearly independent, where 

ao = (50) 

ai = tk-2 - tk-3 = Tk-2 (51) 

0-2 = tk-i - tk-3 = Tk-i + Tk-2 (52) 



J, = 



-b 
a 
A 



(6>0) 



(53) 



with zq,T defined as before. 

The problem is treated in three-dimensional space, and the geometric struc- 
ture can be described as follows. The generic vector Y{a) is written as 



Y{a) — [exp{aa) cos{ba), exp{aa) sin{ba) , exp{Xa)]' 
It is therefore, a spiral on the surface of revolution 



(54) 



(55) 



whose form is determined by the real part of the system eigenvalues. 

The vectors Yq,Yi,Y2 have their origin at the point (0,0,0) and their ends 
at the points ^(ao)! ^^(^1), ^^(^2) of the parametric curve. 

Heuristically we can imagine the same vectors as before (n — 2) pointing 
upwards from the {X, Y) plane as they have a third component on the Z axis. 

From this geometric interpretation we can study some interesting cases of 
linear dependence. 



(a) If the sampling-period sequence is chosen in such a way that 



6a,; 



(i = 0,1,2) 



(56) 



according to Fig. 1, the vectors Yo,Yi,Y2 will be coplanar (in a plane 
containing the Z axis) and therefore they will be linearly dependent since 
every vector can be written as a linear combination of the other two. 
In physical terms, this can be regarded as a resonance of the sampling 
sequence with a pair of complex eigenvalues. This situation can be gener- 
alized to higher-order models if the equation (|56|) holds for the imaginary 
part of any pair of complex eigenvalues. 

(&) If the real eigenvalue and the real part of the complex pair are equal then 
the revolution surface is a cone and there will be linear dependence if at 
least two vectors are on the same generator (see Fig. 2). Analytically, 



b{aj — ai) — 2Tr 



(57) 
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for any i, j such that < i < j < 2. 

Physicahy, this situation can be regarded as a resonance of ah the eigenval- 
ues with the time interval between two sampling instants, not necessarily 
consecutive. 

This situation can be generalized to higher-order models if the eigenvalues 
are such that 

(i) the real parts are the same, and 

(ii) all the imaginary parts bj {j = I, . . . ,q) satisfy ((57)) . 

(c) This is the general case where there is linear dependence because the vec- 
tors Yi are contained in an arbitrary plane, which obliquely intersects the 
surface of revolution passing through the origin. 

The intersection of the oblique plane and the surface of revolution is a closed 
curve r. Therefore the spiral will intersect such a curve in only two points 
for each rotation of 27r radians. From this interpretation, different geometric 
structures can be obtained, depending on the parameter values. In fact, we have 
the following possibilities: 



Now A > (A < 0) implies that the spiral goes up (down) as far as a — oo; 
and a > (a < 0) implies that the surface of revolution expands (shrinks) when 
we move up the Z axis. 

Case 1: A > 0, a > 

From examination of Pq and Pi the normalized projections of Yq and Yi on 
the (X, Y) plane in Fig. 4, it is clear that if Y2 is contained in the oblique plane 
then its normalized projection on the (X, Y) plane is situated on the arcPoPi = 
bai (here Pq and Pi denote the ends of the vectors Po and Pi respectively). 
Therefore a sufficient condition that guarantees the linear independence of the 
three vectors is that the rotation ba2 is not part of such an arc. 

In this way, PqPi and the corresponding multiples of 2tt are the forbidden 
intervals, and consequently PiPq and the corresponding multiples of 2tt are the 
allowed ones. 

It has to be pointed out that inside the forbidden intervals there would 
only be two intersection points of the spiral with the curve F for each rotation 
through angle 2tt. However, we reject the whole interval since we intend to 
find sampling intervals in which the linear independence of the vectors Yi is 
automatically guaranteed without analytical computations. 

The most favourable case corresponds to a rotation PqPi as short as possible, 
because this would increase the length of the allowed interval. Reciprocally, the 
most disfavourable case corresponds to a rotation PqPi, with angle close to n. 



A > 




A < 
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For values of a such that 



\Y{a)\ > \Y{a*)\ 



(58) 



where Y{a*) is the vector of maxnTium modulus with its origin at the point 
(0, 0, 0) and its end on the curve F, the choice of the third vector is completely 
arbitrary, since Y2 can never be on the oblique plane. Therefore the number of 
forbidden intervals is finite and depends on the previous relation. 

For the case A > 0, a < the allowed intervals are the same as before. 

Case 2: A < 0, a < 

We proceed in the same way as before. Figure 5 shows the projections of 
Yq,Yi and the curve F on the {X,Y) plane. The diameter i?ii?o is parallel to 
the chord joining the end of the normalized projections Pq and Pi. For each 
rotation through angle 2-k, the spiral will intersect the curve F in only one point 
at PiRi, and the same will happen at RqPo- Therefore a sufficient condition to 
guarantee the linear independence of the three vectors in that the rotation 60:2 
is not part of such an arc. Consequently, PqPi and RiRq and the corresponding 
multiples of 27r radians would be the allowed intervals. 

It should be pointed out that the arc is allowed only for the first 

rotation. The most favourable case corresponds to a rotation PqPi = bai, 
with angle as close to tt as possible. Conversely, the most unfavourable case 
corresponds to a rotation with angle as small as possible. 

In this situation the number of forbidden intervals is not finite, since there 
are always two intersection points for each rotation of angle 27r. 

For the case A < 0, a > the allowed intervals are the same as before. 

Comparing these results with those of (9j, we can see that for the latter the 
allowed interval, in rotational form, can be expressed as 



which, according to the geometric interpretation, is clearly a shorter interval. 
In both cases the conditions are sufficient. 

The only interesting cases for second- and third-order models are those con- 
sidered in §§ 4.1 and 4.2. 

For systems without complex eigenvalues, the linear independence of the 
vectors Yi is automatically guaranteed for any arbitrary choice of the sampling 
instants [S]. 

4.3 General considerations for higher-order models 

From the fourth-order upwards, the situation becomes more and more complex 
because we do not have the geometrical insight given by the plane or three- 
dimensional space. However, from the basic structures developed for and M^, 
we can make some general considerations that simplify the choice of sampling 
instants in higher-order models. 



ba2 < TT 



(59) 
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In order to clarify these ideas, we are going to consider a 4th-order model 
(n = 4) with one pair of complex eigenvalues and one pair of real eigenvalues: 

Ai, A2 e M, a + jbeC (60) 

The parametric curve is 

Y{a) ~ [exp{aa) cos{ba),exp{aa) sin{ba),exp{Xia),exp{X2a)]' (61) 
which can be decomposed into 

ci{a) = [exp{aa) cos{ba),exp{aa) sin{ba),exp{Xia)] (62) 

C2(a) = [exp{aa) cos{ba),exp{aa) sin{ba),exp{X2a)] (63) 

and we can study the evolution of each curve separately in the same way as 
before. 

The normalized projections of the vectors 

(ci(ai)), (c2(ai)) (2 = 0, ...,3) 

on the {X, Y) plane will coincide, since the only element varying is the third 
component and not the angle of rotation. 

Take {ao, ai, 02, 03) such that each subset of three elements in either of the 
sets (ci(ai)) (c2(ai)) is linearly independent. Then we have that each subset of 
three elements in 1^ (z = 0, . . . , 3) is linearly independent. Therefore it suffices 
to check the linear dependence of any one vector with respect to the other three. 

In analytical terms, this means that 

Miciiaj) ^ M2C2{aj) (64) 
where Mi and M2 are 3x3 matrices of general form 

If ([M)) does not hold then the linear independence of the vectors Yi {i = 
0, . . . , 3), obtained through the manipulation of 3-dimensional vectors, is guar- 
anteed. Therefore we have been able to reduce the problem dimension by one. 

The method can be generalized to arbitrary order n, although the problem 
becomes more and more complicated as the order of the model is increased. 

If the multiplicity of the poles is greater than 1, the third component of the 
vectors Ci{a) is a polynomial expression with trigonometric functions, which will 
give a more complicated curve. 

If there is more than one pair of complex eigenvalues, that with the great- 
est imaginary part (maximum eigenfrequency) will be the main pair; since it 
determines the rotation with greater angle. 

When the dimension of the model is very large, it is more convenient to use 
the shorter intervals found in the literature or to reduce the order of the model 
and then apply the techniques developed to the reduced model. 
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5 Particularization to the periodic case 



The aperiodic model developed in previous sections is also valid for the periodic 
case. In this situation, the general expressions and the geometric interpretation 
can be simplified. 



5.1 Periodic model 

Once the sampling period Tq has been fixed, the multivariable functions fi,gj 
are reduced to constant coefficients , bj . 

It can be demonstrated by induction that for the one- variable case the coef- 
ficients ai can be written as 

a, = (-ir-iE0,i...</),, (66) 
where the sum is taken for 1 < jp < n and 1 < p < i {i = 1, Idots, n) 

(I'jtiTo) = exp{XiTo) (67) 

{Xi is an eigenvalue of the continuous system). If the multiplicity of Ai is 
mi > 1 then there are rrii functions 4>i which are all equal. The a.^ correspond to 
coefficients with opposite sign of the polynomial A{z) in the Z-transfer function 

C(>, ^ |il (68) 

Indeed, 

A(z) = l + a'iZ-i + ... + a;,z"" (69) 

with 

a'i = -ai (70) 

It is well known that between the coefficients and the roots of a polynomial 
there is a variational relation like that given by (|66|) . 
According to p4l) . the coefficients bj are 

n 

bk-j = -'^a^h{tk-i-tj) (j = k,k-l,...,k-n + l) (71) 

1=0 

with oq — — 1, fc > n. Equations (I66p and (|7ip represent an alternative 
approach to that presented by the Z-transform in order to compute discrete- 
model coefficients from the weighting function. 

The formulation corresponding to the case with zeroth-order hold can be 
also transcribed to the periodic version by means of 
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Table 1: Values of the parameters for different sampling periods 



To 


2 


4 


6 


8 


10 


12 


hi 


0.0026 


0.0186 


0.0510 


0.0989 


0.1586 


0.2260 


b2 


0.0092 


0.0486 


0.1086 


0.1718 


0.2257 


0.2643 


h 


0.0018 


0.0078 


0.0139 


0.0174 


0.0181 


0.0167 


ai 


2.2549 


1.7063 


1.2993 


0.9953 


0.7668 


0.5938 




-1.689 


-0.958 


-0.547 


-0.314 


-0.182 


-0.106 


as 


0.4203 


0.1767 


0.0742 


0.0312 


0.0131 


0.0055 


E&i = 1 + Ea.; 


0.0139 


0.0750 


0.1736 


0.2882 


0.4025 


0.5071 



5.2 Influence of sampling period on the discrete model 
parameters 

To discuss the effect of sampling period on the absolute values of parameters 
we give the following example. 

We consider the continuous transfer function 

K 

^^'^ ^ {l+Tis){l+T,s){l + ns) ^'^^^ 
with zeroth-order hold. The external representation is 

3 3 

Vk = ^ atVk-i + ^ b'^Uk-j (73) 

i=l j=0 

where and bj can be obtained according to (l66l) and ([7T|) with 

UTo) = expi-^To) {i = 1, 2, 3) (74) 

For the following values of the continuous model parameters 
K ^1, Ti = 10s, T2 = 7.5s, T3 = 5s 

the values of the discrete-model coefficients are shown in Table ([1]) for dif- 
ferent sampling periods Tq. 

These results agree with those obtained by Isermann |5J. Because of the 
sign of Ti and the general form of , bj , the magnitudes of the parameters 
decrease (in absolute value) and those of bj increase with increasing sampling 
period Tq- 

If Ti had the opposite sign the parameter behaviour would be entirely dif- 
ferent. 

For a small sampling period Tq = Is 

bi < la.l, Sbi < laJ 
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This is why small errors in the estiraated parameters can have a significant 
influence on the input-output behaviour of the model. Indeed, T,bi depends on 
the 4th or 5th place of bi after the decimal point. 

If the sampling period is chosen too small, ill-conditioned matrices result, 
which leads to numerical problems. On the other hand, if the sampling period 
is chosen too large then the dynamical behaviour is described inexactly. For 
Tq = 10 s the model is practically reduced to second order because 

as < 1 + Elflil, 63 < 

and for even greater sampling periods we get a first-order model. 

A proper choice of sampling interval in most cases is not critical, because 
the range between too-small and too- large values is relatively wide [5] . 

5.3 Geometric particularization 

The geometric interpretation of §3 can be particularized to the periodic case. 
We are going to consider again the third order model because it is the most 
significant. For the cases previously presented, we get the following results. 

(a) The vectors Yi are coplanar if and only if 

bTo = IT (75) 

(6) This is included in the previous case. 

(c) This situation never occurs, since different rotations of the vectors Yi would 
imply different sampling periods, which is not possible in the periodic case. 

We can see that the transition from the aperiodic case to the periodic case 
is a simple particularization. In the opposite sense, the procedure is much 
more complicated, because there are many aperiodic situations that have no 
equivalent in the periodic version. 

In analytical terms, the coefficients a^, and consequently the bj, are perfectly 
defined for all sampling periods To, according to ([55]) and ([7T|) . We recall that 
the Z-transform imposes no constraints on the choice of Tq. Conversely, in 
the aperiodic formulation there are some sampling sequences that are forbidden 
for the model developed. This confirms the complexity of the aperiodic case 
compared with the periodic one. 

5.4 CoefRcients of the discrete model for systems with 
dead time 

We are going to consider the same system with zeroth-ordcr hold as before, but 
also with dead time Td- The convolution expression that reflects this situation 
is 
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fc 

yk = ^ h{tk -ti~ Td) Xi (76) 

i=0 

with Xi and h{t) defined as before. For 

{p - l)To <Td< pTo (77) 
with p ^ the general input-output model can be written as 

71 n 

Vk = ^ aiUk-i + ^ bjUk^j-p (78) 

i=l j=0 

In order to determine the coefficients and 6j, the procedure is similar to 
that previously developed. Indeed, now we have the expressions 

n 

J2 a» [ Htk-^ - tj - Td) - h{tk-., - tj+i -Td)]+ bk-j = (79) 
with (j = k — p,k — p — — p — n). 

n 

J2a^[Htk-^-tJ-Td)-h{tk-^-tJ+l-Td)]^0{J ^k-p-n-l,...,0) {80) 

4=0 

and we see that for j = fc — — p+1 

h{tk - -Td) = (81) 
There are two different cases. 

(i) Td — pTq, which implies that the coefficients are not modified by the 

dead time; the coefficients bj are not modified either, except that these 
coefficients multiply the impulse inputs Uk-j-p {j = 0, 1, . . . ,n), rather 
than the inputs Uk-j {j = 0, 1, . . . , n) for the case without dead time. 

(ii) (p— l)To < Td < pT(), which implies that the coefficients are not modified 

by the dead time; the coefficients bj are now functions of Td whose general 
expression is given by (I79p and which multiply the inputs Uk-j-p (j = 
0,1,. ..,n). 

6 Conclusions 

A general aperiodic model for linear time-invariant SISO systems has been de- 
veloped. The model also covers more general cases such as systems with zeroth- 
order hold and dead time. The formulation considered stresses the importance 
of the sampling period against other system parameters. In this way, such sys- 
tems have an additional element for analysis and manipulation. The periodic- 
sampling case appears as a simple particularization of the general procedure. 



18 



The results obtained are simplified and the use of tables of Z-transfornis are 
avoided. For every sampled system, in a periodic or aperiodic way, there will 
always be sampling-period sequences more or less adequate according to the 
general characteristics of the process under study. In the aperiodic case, there 
will also be some restrictive conditions on these sequences, although it has also 
been possible to give strategies for some special cases. 
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